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Abstract 

In this work, we recalculate the charmless pure annihilation decays Bg — )• vr+vr^ and — )■ 
Clf K+K- by using the perturbative QCD (pQCD) factorization approach, and compare the pQCD 

predictions with currently available experimental measurements. By numerical calculations and 
phenomenological analysis we found the following results: (a) one can provide a consistent pQCD 
interpretation for both the measured Br{B^ — t- tt^tt^) and Br(B^ — t- K~^K~) simultaneously; 
'■^^ (b) the pQCD predictions for Br{B^ — t- tt^tt^) obtained by different authors are well consistent 

with each other; (c) our new pQCD prediction for Br{B^ — )• K^K~) agree well with the 



measured values from CDF and LHCb Collaboration; and (d) the CP-violating asymmetry 
-^Cp{B^d K+K-) ^ 19%, which is large and may be detected at the LHCb and future 
Super-B factory experiments. 
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I. INTRODUCTION 



Among the two-body hadronic B meson decays, the pure annihilation decay modes, 
such as — ?► TT+TT^ and — t- K^K^ decays, are specific in several respects. They can 
occur only through the annihilation diagrams in the standard model (SM) because none 
of the quarks(anti-quarks) in the final states are the same as those of the initial B meson. 
And consequently, they are rare decay modes with a branching ratio at the level of 10~^ 
or less as generally expected. Such decays play very important role in understanding the 
annihilation mechanism and determining the strength of the annihilation contribution in 
B meson charmless hadronic decays, and therefore have been studied intensively by many 
authors [l-lol| in spite of the great difficulties in both the theoretical calculation and the 



experimental measurements. 

In the experiment side, both Br{B^ — t- tt+tt^) and Br{B^ — t- K^K^) are measured 
very recently due to their rareness. At the spring and summer conference of 2011, CDF 
111 and LHCb fl2t| collaboration reported their first measurement of the decay rates 

rwrO_^^+ -^_ / (5.7±1.5(stat.)±1.0(.yst.))xl0-^ CDF [11] 
vr TT J - I (9 8+2.3(^^^^ ) ^ l.l{syst.)) x 10^^ LHCb [12], 

RWR0^/<'+/<'-l-/(2-3±l-0(stat)±1.0(sy.t.))xl0-^ CDF [11], 
K J-|(i.3+o.6(^^„^_)±0.7(si/sl))xlO-^ LHCb [12]. 

The statistical significance of LHCb measurement reaches 5.3cr for Bg — ?► tt+tt^ decay, 
which means a observation for the first time. 

In the theory side, we know that it is very hard to make a reliable calculation for pure 
annihilation decays of B mesons. In the QCD factorization (QCDF) approach (l3|. for 
example, one can not perform a real calculation for the annihilation diagrams due to the 
end-point singularity, but have to make an rough estimation by parameterizing the annihi- 
lation contribution through the treatment dx/x ^ Xa = (1 + Pa^^'^) In ^ l^l, Isl, or by 
using an effective gluon propagator — )■ 1/ (A;^ + Mglk"^)^ to avoid enhancements in 
the soft endpoint region [7]. Of course, such parameterization will produce large theoret- 
ical uncertainties. For 5° — )■ tt+tt^ and B^ — )■ K^K~ decays, the theoretical predictions 
based on the QCDF approach as given for example in Refs.jsHS] are the following: 

( 0.24 X 10-^ [5,6] 
5r(5° ^ 7r+7r~) = <^ (1.24 ± 0.28) x 10-^ [7] (3) 
( (2.6 ±1.0) X 10^^ [9] 



Obviously, the QCDF predictions in Refs. [5|-|8| are much smaller than the measured results 
for Br{B2 — )■ tt+tt"), while smaller or close to the measured ones for Br{B^ — )■ K^K~) 
in Ref. j5[ and Ref. [Sj, respectively. 



After CDF's report of the evidence of Bg tt+tt^ decay, the author of Ref. [10 
reinvestigated the role of annihilation topology in the QCDF approach and found that 

(1) the CDF measurement of Br{Bs tt+tt^) implies a large annihilation scenario with 
Pa around 2 instead of ~ 1 preferred by all previous studies in QCDF approach jHI,^, 13]; 

(2) if one assumes universal annihilation parameters pa and (pA for all B^ ,, — )■ PP decay 
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modes, one can not provide predictions being consistent with all well measured decay£]; (3) 
one possible way to solve this problem is to use different (p^, (^a) for different decays, which 
however means that the predictive power of QCDF approach becomes rather limited. In 
short the studies in Ref. [lo| tell us that it is very hard to give a consistent QCDF 
interpretation for Br{Bs — > vr+vr") and other well measured Bd^s — )■ PP decay modes 
simultaneously. 

In the perturbative QCD (pQCD) factorization approach 17|. however, the situa- 



tion becomes rather different. Here, the pure annihilation decays of B/Bg meson can be 
calculated perturbatively by employing the Sudakov factors to smear and then to strongly 
suppress the end-point singularity. In the pQCD factorization approach, for example, the 
endpoint divergence of the factorizable emission diagram Fig. 1(a) and 1(b) in Ref. [isj 
are regulated by introducing the transverse momentum kf^, i.e. 



1 1 1 



F pI-MI M|xiX3(l-X3) 

1 



(5) 



where k"^ = {ki — k^y and pi = (Pi — k^,)^ is the momentum of the gluon propagator 
and 6-quark propagator respectively. It is easy too see that the end-point divergence for 
Xi = and x^ = (0, 1) are removed effectively by introducing small but non-zero k?^. 

For — >■ 7i~^7i~ decay, it was calculated by employing the pQCD factorization ap- 
proach in 2004 and 2007 Q, respectively. In Ref. 0, we obtained the first pQCD 
prediction for the decay rate: 

fir(fi° TT+TT") = (4.2 ± 0.6) X 10~^. (6) 

In 2007, Ali et al, j3| made a systematic calculation for all B^ — )■ PP, PV, VV decays in 
the pQCD factorization approach and found that 

5r(5° ^ TT+TT-) = (5.71};^) x 10"^. (7) 

These two pQCD predictions at leading order (LO) are well consistent within la error and 
confirmed by CDF and LHCb measurements as shown in Eqs. f lTP]) . The small difference 
for the predicted decay rates between Ref. ^ and jsj comes from the fact that a little 
different input parameters and distribution amplitudes(DA's) of vr and Bs meson were 
used in two studies. 

In Ref. 0, by employing the pQCD factorization approach, we studied the Bg — )• PP 
decays with the inclusion of partial next-to-leading order (NLO) contributions, coming 
from the QCD vertex corrections, the quark-loops, the chromo-magnetic penguins and the 
usage of the NLO Wilson coefficients instead of the LO ones. For the pure annihilation 
decay Bg — > tt+tt^, it does not receive the NLO contributions from the QCD vertex 
corrections, the quark-loops and the chromo-magnetic penguins. The leading order pQCD 
prediction is Br{Bg — )■ tt+tt" -) = (7 ±2.5) X 10"^ while it becomes (5. 712:2) x 10"^ when 
the NLO Wilson coefficients Ci{Mw), the NLO renormalization group evolution matrix 



^ The corresponding QCDF predictions for Br{Bs K^K ) and Br{Bd K'^K^) are twice larger 
than the experimental measurements [lo| . 
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U{t, m, a) [191 and the as(t) at two- loop level were employed in the numerical calculation 
i- 

For — )■ K^K decay, the known pQCD prediction for its branching ratio was given 
in 2001 [ij 

Br{B^^ K+K-) = 3.27 x 10~®; 5r(5^ -> K+R-) = 5.90 x 10~^ (8) 

which is much smaller than the measured value as given in Eq. ([2]) by roughly a factor of 
three, in other words, a large discrepancy between the data and the theoretical prediction 
based on the pQCD factorization approach for B^ — )■ K^K~ decay. 

It is necessary and interesting to check if one can provide a consistent pQCD interpre- 
tation for both the measured Br{B^ — )■ vr+vr") and Br{B^ — )■ K^K~) simultaneously? 
In this paper, by employing the pQCD factorization approach, we recalculate the pure 
annihilation decays B^ — ?► 7r+7r~ and B^ — )■ K^K~ with the usage of the same set of input 
parameters and wave functions for the mesons involved, in order to check if the new data 
from CDFF and LHCb can be understood in the pQCD approach. Our studies will be 
helpful to determine the strength of penguin-annihilation amplitudes (2l| . 

The paper is organized as follows. In Sec. HIl we give a brief review about the theoretical 
framework of the pQCD factorization approach and the wave functions for B^/B^ and 
TT, K mesons involved. We perform the perturbative calculations for considered decay 
channels in Sec. IIIH while the numerical results and phenomenological analysis are given 
m Sec. HVl A short summary also be given in Sec. IIVI 



II. THEORETICAL FRAMEWORK 

In the pQCD approach, the decay amplitude A{Bg — M2M3) with q = {d, s) can be 
written conceptually as the convolution, 

A{Bg ^ M2M3) ~ jd''hd%d% Tr[CmB{ki)<^M,{k2)^M,{k3)H{h,h,h,t)] , (9) 

where fcj's are momenta of light quarks included in each meson, and "Tr" denotes the 
trace over Dirac and color indices. In the above convolution, C{t) is the Wilson coefficient 
evaluated at scale t, the function H{ki, ^2, ^3, t) describes the four quark operator and the 
spectator quark connected by a hard gluon. The wave function and describe 

the hadronization of the quark and anti-quark in the Bg meson and the final state light 
meson Mj. 

We treat the Bq meson as a heavy-light system, and consider the Bq meson at rest for 
simplicity. Using the light-cone coordinates the Bq meson momentum Pb and the two 
final state meson's momenta P2 and P3 (for M2 and M3 respectively) can be written as 

P« = -|(1,1,0t), P2 = ^(l-r3^r2^0T), = ^{rl\ - rl^^) , (10) 

where rj = nii/MB- For the final state light mesons made up with {u,d,s) and the 
corresponding anti-quarks, the ratio r2 and are small and will be neglected safely. 
Putting the quark momenta in Bq, M2 and M3 meson as ki, ^2, and k^, respectively, we 
can choose 

fci = (a;iP+,0,kiT), fc2 = (x2P2+,0,k2T), ^3 = (0,X3P3-,k3T). (11) 
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Then, the integration over , k2 , and in eq.(|9]) will lead to 



A{B^ ^ M2M3 



dxidx2dxsbidbib2db2bsdb3 



•Tr 6i)$Aft(x2, &2)$Af3(x3, bs)H{x„ h, t)St{xi) e-^W] 



(12) 



where 6j is the conjugate space coordinate of ki^. The large double logarithms (In^ Xi) on 
the longitudinal direction are summed by the threshold resummation, and they lead to 
the first Sudakov factor St{xi) which smears the end-point singularities on Xi [l7|. The 
Sudakov resummations of large logarithmic corrections, such as the terms proportional to 
as log^[Q/kiT] ( Q ~ m^), to the Bg and two final state meson wave functions will lead to 
the second Sudakov factor e"'^*-*-' = e~^'^^^^ .^-^M^ii) .^-SM^{t)_ xhese two kinds of Sudakov 
factors can together suppress the soft dynamics effectively (l7l |. 

In the momentum space, the light-cone wave function of Bq meson can be defined 



as 



'+"^ijj750i?„(^) 



(13) 



where P is the momentum of the Bg meson, k is the momentum carried by the light quark 
in Bq meson, and (pBq is the corresponding distribution amplitude. 

For the B/Bg mesons, the distribution amplitudes (psix, b) in the b space can be written 



as 



14-16 



and 



0b(x, b) = Nbx'^{1 — xY exp 



(j)B,{x,b) = Nb^x'^{1 - x^exp 



xiriB 



ulb^ 



(14) 



1 / xniB, 

2 V ^B. 



Bs 



(15) 



where the normalization factors Nb^^-^ are related to the decay constants /s,^, through 



dx(t)B(,){x,b = 0) 



Ib, 



2^6 



(16) 



Here the shape parameter ujb has been fixed at 0.40 GeV by using the rich experimental 
data on the B mesons with Jb = 0.19 GeV. Correspondingly, the normalization constant 
Nb is 91.745. For Bs meson, considering a small SU(3) symmetry breaking, since s quark 
is heavier than the u or d quark, the momentum fraction of s quark should be a little 
larger than that of m or d quark in the B mesons, we therefore adopt the shape parameter 
ubs = 0.50 GeV ^] with /^^ = 0.23 GeV, then the corresponding normalization constant 
is A''^^ = 63.67. In order to analyze the uncertainties of theoretical predictions induced by 
the inputs, we can vary the shape parameters Ub and ub^ by 10%, i.e., Ub = 0.40±0.04 GeV 
and ubs = 0.50 ± 0.05 GeV, respectively. 

For the n"^ and K"^ mesons, we adopt the same set of distribution amplitudes 4>^K{xi) 

IM)- 



and 



Xj 1 clS defined in Refs. 



3/7r(j^) 



x(l — x) 



l + a^(^)cf (t) + a2(^) 



(17) 
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Bs/B 






id) 



FIG. 1. The typical annihilation Feynman diagrams for B'^ — )• ir^ir and — )• K^K de- 
cays, (a) and (b) are factorizable diagrams; while (c) and (d) are the non-factorizable ones. 



\{K) 



X] 



f-K(K) 



(18) 
(19) 



where t = 2x — 1, Pt,{k) = fn-KiK) /n^'o^^^ are the mass ratios ( here mg = m\/{mu + md) 
and = m|-/(ms + mrf) are the chiral mass of pion and kaon), a^'^ are the Gegenbauer 
moments, while C^(t) are the Gegenbauer polynomials 

Cf/^(t) =3t, 



(3-30t2 + 35t^) , Cf'(t) = ^(l-14t2 + 21t^) 



15 



(20) 



Under the replacement oi x ^ 1 — x, only Cf''^(t) will change its sign, others remain 
unchanged. 



III. PERTURBATIVE CALCULATION IN THE PQCD APPROACH 

In the pQCD factorization approach, the four annihilation Feynman diagrams for 
Bs — )■ vr^vr" and B^ — )• K^K~ decays are shown in Fig.l, where (a) and (b) are fac- 
torizable diagrams, while (c) and (d) are the non-factorizable ones. The initial h and s{d) 
quarks annihilate into u and u pair, and then form a pair of light mesons by hadronizing 
with another pair of dd {ss) produced perturbatively through the one-gluon exchange 
mechanism. Besides the short- distance contributions based on one-gluon-exchange, the 
qq pair can also be produced through strong interaction in non-perturbative regime (final 
state interaction(FSI), for example). FSI effects in considered decays have been assumed 
rather small, we do not consider them here. 

We will adopt (F^^, F^^, F^^) and (M^^, M^^, M^^) to stand for the contributions of 
the factorizable (Fig. 1(a) and 1(b)) and non-factorizable (Fig. 1(c) and 1(d)) annihilation 
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diagrams from the {V - A){V - A), {V - A){V + A) and {S - P){S + P) operators, re- 
spectively. By making the analytic calculations we obtain the following decay amplitudes 
for both B° tt+tt" and S° K+K~ decays: 

From the factorizable annihilation diagrams Fig. 1(a) and Fig. 1(b) we have 

(i) {V - A){V - A) operators: 

P 1 POO 

F^^ = IQttCfMI^ / dx2dx3 / b2db2hdb3 
Jo Jo 

X { [X20^(X2)0^(X3) + 2r2r3 (0f (Xs) - 0^(X2)) 0f (Xg) 
+2r2r3X2 {4>2{^2) + 0^(3^2)) 0f (3:3)] ■ ha{x2, X3, 62, ^s) -Ea(ta) 

+ [{X3 - l)0^(x2)(/)3 (X3) - 4r2r30f (a;2)0f (X3) 

+2r2r3X30^(x2) (0f (xg) - ^n^^s))] ■ h{x2,X3,b2M) Eaitb)] , (21) 

(ii) {V - A){V + A) operators: 

pLR ^ pLL ^ (^22) 

(iii) (5 — P){S + P) operators: 

F^P = ?,2txCfMI^ j dx2dx3 / b2db2b3db3 
' Jo Jo 

X { [2r30^(a;2)0f (X3) + r2X2 (0f (^2) - 0^(x2)) 0^(0:3)] 

■ha{x2,X3,b2,b3) ■ Ea{ta) 

[2r20r(x2)0^(x3) + (1 - X3)r3</)^(X2) (0f (X3) + ^n^^s))] 
■hb{x2,X3,b2M) ■ Ea{tb)} , (23) 



where r2 = 'm2/mB^, T3 = m^/niBq with g = (c/, s) for 5^ or 5° decays, and Cj? = 4/3 
is a color factor. The explicit expressions for the convolution functions Egjtn^h), the hard 



scales ta,b, and the hard functions ha^bi^i, bi) can be found for example in Ref. [18|, |20 
From the non-factorizable annihilation diagrams Fig. 1(c) and Fig. 1(d) we have 

64 

M^^ = —^TrCpM"^ / dxidx2dx3 I bidbib2db2 4>B^{xi^bi) 
V 6 'Jo Jo 

X { [(1 - X3)ct>^{x2)cl>i{x3) + (1 - a;3)r2r3 (0^ (0:2) + 0^X2)) (^g^'la^s) - 0r(a;3)) 

+X2r2r3 (0f (X2) - 02 (a;2)) (0f (a^s) + 03 (a^s))] ■ hc{xi,X2, X3, bi, 62) Ec(tc) 

+ [-a;20^(a;2)0^(a;3) - 4r2r30f (x2)0f (X3) 

+ (1 - X2)r2r3 (0f (X2) + 0^(X2)) (0f (X3) - (Plixs)) 

+X3r2r3 (0f (X2) - 02 (a;2)) (0f (a;3) + 03 (a;3))] 

■hd{xi, X2, X3, 61, 62) Ec{td)} , (24) 
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M 



LR 



64 

—j=TiCpM'^ I dxidx2dx3 / bidbib2db2 4>Bq{xi,bi) 
v6 "Jo Jo 

X { [x2r2 {^2(^2) + 4'^ix2)) (Ptix-i) - (1 - a;3)r30^(x2) (0f (xg) - 03 (a^s))] 

■hciXi,X2,X3,bi,b2) Ecitc) 

+ [(2 - X2)r2 (0f (X2) + 0^0:2)) 03 (a^s) - (1 + a;3)r30^(x2) (0f (X3) - 03 (a^s))] 
■hd{xi,X2,X3,bi,b2) Ec{td)} , (25) 



M' 



SP 



bidbib2db2 <PbAxiM 



r-l poo 

dxidx2dx3 I 
'0 Jo 

X { [x20^(a;2)0^(a;3) + a;2r2r3 (0^X2) + 0^(2:2)) (0^(3:3) 

+ (1 - X3)r2r3 (0^(X2) - 0^(X2)) (0f (X3) + 0^X3))] 

■hc{xi,X2,X3,biM) Ec{tc) 
+ [-(1 - X3)0^(X2)0^(X3) -4r2r30f(a;2)0f (0:3) 
+a;3r2r3 (0^(3:2) + 0^(x2)) (0^(0:3) - 0r(x3)) 
+ (1 - X2)r2r3 (0^X2) - 0^(X2)) (0f (X3) + 4>1{X3))] 

■hdixi, X2, X3, 61, 62) E^{td)} , 



01(^3)) 



(26) 



where r2,3 and Cp are defined in the same way as in Eqs. (12111231) . Again, the explicit 
expressions of the functions -Ec(^c,d) and hc^d, and the hard scales tc,d can be found in 
Ref. [l8l,[23. 

Because of the isospin symmetry, the contributions to both — t- tt+tt^ and — t- 
K^K' decays from the factorizable annihilation diagrams Fig. 1(a) and Fig. 1(b) cancel 
each other. The total decay amplitudes for the considered decays are therefore written 



as: 



A{Bl ^ TT + TT-) = V:,Vus C2M 



LL 



M 



LL 



M 



SP 



(27) 



^(S° ^ K^K-) = V:,V^d C2M 



LL 



C4 + — -Cg + Clo 



C4 + Cq + Cg — 2^10 



M 



SP 



■ (28) 



The expression of decay amplitude in Eq. fl27|) is equivalent with those as given in Refs. j2|, 
a by a proper transformation between M^^ and M^^. 



IV. NUMERICAL RESULTS AND DISCUSSIONS 

Now it is straightforward to calculate the CP-averaged branching ratios and CP- 
violating asymmetries for the two considered decays. In numerical calculations, central 
values of the input parameters will be used implicitly unless otherwise stated. The QCD 
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scale (GeV), masses (GeV), decay constants (GeV), and Bq meson lifetime (ps) being 
used are the following [24] 

Aqcd = 0.25 , mw = 80.41 , m^o = 5.2795 , Mb, = 5.3663 ; 
= 0.14, TTiK = 0.494, = 0.13, fx = 0.16, 
tbo = 1.525 ps, tb, = 1.472 ps. (29) 



As for the CKM matrix elements, we use [24 



A = 0.2253 ±0.0007, A = O.SOSto^oll P = 0.132l°:°22, = 0.341 ± 0.013. (30) 

For the Gegenbauer moments and other relevant input parameters, based on the works 
of (221,123 



a^ = 0, af = 0.06 ±0.03, = 0.35 ± 0.15, af = 0.25 ± 0.10, 

= -0.015, af = 0, p^ = m^/ml, pK = mK/mQ, 
r/3 = 0.015 ± 0.005, W3 = -3.0 ± 1.0 (31) 

with the chiral mass = 1.4 ±0.1 GeV, and mf = 1.9 ± 0.2 GeV. In order to check 
the theoretical errors induced by the uncertainty of the Gegenbauer moments we vary af" , 
as'^, r/3 and U3 in the range of af = 0.06 ± 0.03, = 0.35 ± 0.15, af = 0.25 ± 0.10, 
r/3 = 0.015 ± 0.005 and Us = -3.0 ± 1.0. 

From the decay amplitudes, it is easy to write down the corresponding branching ratio: 

Br{B ^ PP) = TB \A{B ^ PP) f , (32) 

where A{B — > PP) is the decay amplitude as defined in Eqs. (27,28). 

By using the analytic expressions for the complete decay amplitudes and the input 
parameters, we calculate the branching ratios and CP-violating asymmetries for both 
considered decay modes. The numerical results are the following: 

i?r(5° ^ vr+TT-) = (5.10tl|^(a-)iOi^tS;itO:l) x 10"^ (33) 

AcpiB', ^ vr+vr-) = (-2.3;°;°(a-)l°i;°il°:J) %, (34) 

Br{B', ^ K^K-) = (l.56;°:-(af )l°:^n?^^S:J^) X 10-^ (35) 

AcpiB', ^ K^K-) = (I8.9t?:^(af )t^:^l?:n?) %, (36) 

where the first error comes from the theoretical uncertainty of the Gegenbauer moments 
a2 = 0.35 ± 0.15 and af = 0.25 ± 0.10, the small theoretical errors due to the variations 
of af = 0.06 ± 0.03, 773 = 0.015 ± 0.005 and U3 = —3.0 ± 1.0 is shown as the second error, 
the third error includes the uncertainties induced by the parameter Ub = 0.40 ± 0.04 GeV 
and ubs = 0.50 ± 0.05 GeV, as well as the uncertainties of uiq = 1.4 ± 0.1 GeV and 
mf = 1.9 ±0.2 GeV, and the last error comes from the uncertainties of the relevant CKM 
elements. It is easy to see that the uncertainties from a^'^ , Ub and ub, dominate the 
theoretical error. 

For Bg — > TT+TT^ decay, the pQCD prediction for its branching ratio in Eg. (l33l) agree 



very well with the measured results from CDF and LHCb collaboration [ll|, |12[ as shown 
in Eqs. ([T]). This pQCD prediction also agree very well with the previous pQCD predic- 
tions as given in Refs. [2|-y. The analytical results for the decay amplitudes obtained in 
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this paper are consistent with those as given in Refs.[2H^. The small difference in numer- 
ical pQCD predictions comes from the difference of the input parameters being used in 
different works. 

For — > K~^K~ decay, fortunately, the pQCD prediction for its branching ratio in 



Eq. ( 135|1 agrees well with the measured results from CDF and LHCb collaboration Il2 
as shown in Eqs. ([2]). 

It is easy to see that the new pQCD prediction in Eq. (j35l) is much larger than the 
previous pQCD prediction as given in Ref. yy]. In order to find the reason for the large 
difference, we checked the relevant analytical expressions as given in Ref. [l| and found 
that those analytical results are consistent with our results after proper transformation: 
x — )■ 1 — x. The large numerical difference between two pQCD predictions comes from the 
fact that (a) the distribution amplitudes of the kaon meson used by Chen and Li [l[ are 
very different from those used in this paper; and (b) some improved Gegenbauer moments 
as given in Ref. [23[ are used in this paper. 

In Ref. [H, only the axial- vector and pseudo-scalar kaon wave functions (pxix) and 
(f)'x{x) were considered: 

^k{x) = ^x{1-x){1 + 0.51(1 - 2x) + 0.3 [5(1 - 2x)'^ - l] } , (37) 
0'i.(a:) = ^x(l-x), (38) 

In this paper, however, besides the leading twist-2 4>k{^) (i-^- the axial-vector 4>k{x) 
in Ref. [l| ), we also take into account the twist-3 contributions from both 0^ and 0^ 
simultaneously. Based on the analytical expressions as given in Eqs. fll7til9p . one can 
obtain the numerical expressions for 0;^(x),0^(x) and 0^(x): 



X 



x{l-x){l- 0.18(1 - 2x) + 0.375 [5(1 - 2^)^ - l] } , (39) 



(f)P(x) = {1 + 0.282(1 - 6x + 6x^) - 0.012 [3 - 30(2a; - 1)^ + 35(2x - 1)^1 } , (40) 
2v6 

^l(x) = {2x - 1) [1 + 0.55 (1 - lOx + lOa;^)] , (41) 

2v6 

by using the central values of the relevant input parameters 4^, Pk,V3 ^tc, as 

given in Eqs.( |29]l3Tl) . 

For the leading twist-2 axial- vector wave function, the used is in the same 

form as (pxix) being used in Ref. jH. The difference of the coefficients of the second 
and third term comes from the variation of the values of the corresponding Gegenbauer 
moments (af,af): (af,af) = (0.17,0.20) in Ref . [iLwhile (af,af) = (0.06,0.25) in this 



paper, based on recent improvements made in Ref. [23|. The difference of the sign of the 
second term in 0^(a;) is resulted from the different assignment for the momentum fraction 
X in Ref. [ij and in this paper: We here use x to denote the momentum fraction of s/s 
quark in the K"^ meson, instead of the u/u quark as assigned in Ref. [l|. The Gegenbauer 
polynomial C^^^(t) = 3t in Eq. (38) will change its sign under the transformation x — )■ 
1 — X. 

In Ref. [H, the authors took (p'^ix) = ^/i^x(l — x) as the pseudo-scalar kaon wave 
function, which was "determined from the data of the B — )■ Kir decays" by Chen and Li, 
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instead of the ordinary 0^(a;) as derived from the QCD sum rule 22|, l23| and shown in 
Eqs. (39). (for more details of the derivation of 0x(a;), see Sec. IV of Ref. (keumOl-b)). 
The (f)'x{x) in Ref. [3] is just the first and leading term of the twist-2 part 4>x{x) and is 
very different from commonly used 0^. 



In Ref. the term (j)J^{x) was absent. All differences in the relevant wave functions 
being used in Ref. [l[ and in this paper lead to the large difference between the pQCD 
predictions for the branching ratio Br[B^ — )■ K^K~) as presented in Ref. [H and in this 
paper. 

Explicit numerical examinations also show that the leading twist-2 0^ provide the dom- 
inant contribution to the magnitude of the decay amplitudes and consequently branching 
ratio Br{B^ K+R-): 

1. When all three terms 0^'^'"^, or only the leading twist-2 term </);^(x), are taken into 
account, we find numerically 

AiB^ ^ K+K-\ - / (-0-31 - 2-2 I) X 10-^ (0^(x) only); 

^ ^ ~ \ (-0.82 - 3.6 I) X 10-^ (All three terms); ^ > 



2. If only the twist-3 term 0^(x), 0^ or both of them are taken into account, we find 
numerically 

A(B« K-K-) = ( (-»f - ^n-T'- iliil m 

^ ' [ (0.06 - 0.27 1) X 10 ^, (0k(x) only); ^ ' 

Br[B' K-K-) = [ 0-0^ " (45) 
^ ' \ 0.01 X 10 (0k(x) only); ^ ' 



It is straightforward to see from the above numerical results that 

1. The leading twist-2 term 0;^(x) provide the dominant contribution to the decay 
amplitude: A = (—0.31 — 2.2 J) x 10~^ if only (px{x) is taken into account, while 
A = (-0.61 - 0.55 /) X 10-5 ( A = (-0.06 - 0.27 /) x 10^^ ) if only ^^{x) ( 0^) 
is taken into account. For the branching ratio, its size would be 10"^ 10"*^ or 10"^ 
if only the term (f)^{x), (p^i^) 'P'xi^) contribute. 

2. The enhancements due to the constructive interference between the three parts also 
play an important role in producing a large branching ratio Br{B^ — )■ K~^K~). 
One can see that the contributions to the decay amplitude A from the three terms 
interfere constructively, which finally leads to a large branching ratio Br{B^ — 
K^K') = 1.56 X 10"^, partially due to the further magnifying effects since the 
branching ratio is proportional to the module square of the decay amplitude A. 

As for the CP- violating asymmetry for the considered decays, Acp{Bg — t- tt+tt^) is 
very small, only about two percent and therefore hardly to be detected even at the LHCb. 
For B^ — K^K~ decay, however, its Acp is relatively large, around 19%, and may be 
detected at the LHCb experiment or future Super-B factory experiments. 
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In summary, by employing the pQCD factorization approach, we here recalculated the 
branching ratios and CP-violating asymmetries of the pure annihilation decays — )■ 
and B° K+K' with the usage of the wave functions based on the QCD sum 
rule [22I, 2^ and the improved Gegenbauer moments [23j. By numerical calculations and 
phenomenological analysis we found the following results: (a) one can provide a consistent 
pQCD interpretation for both the measured Br{B^ — t- tt+vt") and Br{B^ — > K^K~) 
simultaneously; (b) the pQCD predictions for Br{Bg — )■ tt+tt^) obtained by different 
authors are well consistent with each other within one standard deviation; (c) our new 
pQCD prediction for Br{B^ — )■ K^K^) agrees well with the measured values from CDF 
and LHCb Collaboration; and (d) the CP- violating asymmetry Acp{B^g — > tt+tt^) 
—2.3%, may be too small to be detected even at LHCb experiment; (e) Acp{B^ — )• 
K^K~) PS 19%, which is large and may be detected at the LHCb and future super-B 
factory experiments. 
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